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ABSTRACT 

We work in the locally de Sitter background of an inflating universe and consider a massless, 
minimally coupled scalar with a quartic self-interaction. We use dimensional regularization 
to compute the fully renormalized scalar self-mass-squared at one and two loop order for a 
state which is released in Bunch-Davies vacuum at t = 0. Although the fleld strength and 
coupling constant renormalizations are identical to those of flat space, the geometry induces 
a non-zero mass renormalization. The flnite part also shows a sort of growing mass that 
competes with the classical force in eventually turning ofl this system's super-acceleration. 
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I. INTRODUCTION 

For studying quantum fields which are spectators during inflation it is convenient to 
model the background geometry using a portion of the full de Sitter manifold known as the 
open coordinate patch. If the £)-dimensional cosmological constant is A = {D — 1)H'^, the 
invariant element can be expressed, in either conformal or co-moving coordinates, as follows: 

ds^ ^a^(-drf + dx-dx)^-dt^ + a?dx-dx where 0(77) = --^ = e^* . (1) 
^ ^ Ht] 

One might expect that the inflationary expansion of this spacetime makes quantum effects 
stronger by allowing virtual particles to persist longer than in flat space. Indeed, it is simple 
to see that any sufficiently long wavelength (wavelength > 1/H) virtual particle which is 
massless on the Hubble scale can exist forever [1]. However, one must also consider the 
rate at which virtual particles emerge from the vacuum. Classical conformal invariance 
causes this rate (per unit co- moving time) to fall like 1/a, so any long wave length virtual 
particles which emerge become real, but very few emerge [1]. To get enhanced quantum 
effects during inflation therefore requires quanta which are effectively massless and also not 
classically conformally invariant. 

A fascinating example of enhanced quantum effects occurs when a massless minimally 
coupled scalar is endowed with a quartic self-interaction, 

1 A 
C = —-d^(pdy(pg^''\/^ — -^'^'^y/—g + counterterms . (2) 

If this system is released at t = (7/ = —1/H) in free Bunch-Davies vacuum the fully 
renormalized energy density and pressure are [2, 3], 

A XHUI... 2 _3 1~ n + 2 _ ^ 
= 8^+2^(2^^ («) + 9« ^|+0(A^), (3) 

= -8^-2^12^^ ^«) + 3^^^«) + 6£(;^« ^|+0(A^)- (4) 
Hence their sum violates the weak energy condition on cosmological scales, 

+P™ - 4ln(a) + - 1| ^^-j + O(A^) . (5) 

The physics behind this surprising result is that inflationary particle production causes 
the scalar to undergo a random walk such that its average distance from the minimum of 
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the </7^ potential increases without bound in the free theory [4-6], 



= ^Divergent Constant j + ^-^ln(a) . (6) 



After the ultraviolet divergence is removed by renormalization the effect is to drive the scalar 
up its potential, which increases the vacuum energy and leads to (3) and (4). The process is 
self-hmiting for two reasons. First, as the scalar rises up its potential the classical c/?^ force 
pushes it back down. Second, the </7^ curvature associated with being away from (/? = acts 
like a positive mass-squared to reduce the inflationary particle production responsible for 
pushing the scalar away from (p = 

This model is of more than passing interest because the supernova acceleration data are 
consistent with an evolving dark energy equation of state whose current value is less than 
minus one [7-18] . Such a thing is difficult to explain with classical models on account of the 
problem with stability [19-83]. Therefore one is lead naturally to a self-limiting quantum 
effect in a classically stable theory such as (2). 

Although this model gives a value for w + 1 which is far too small to explain the obser- 
vations [3], it does serve to establish that quantum effects can induce a temporary phase of 
super-acceleration. Further, it does so in a wonderfully simple setting for which every issue 
can be settled by exphcit computation. To check stability one flrst computes the self-mass- 
squared M^(x; x') and then studies how it modifies the effective field equations at linearized 
order, 

U^{x) - j d^x'M\x; x')f{x') = . (7) 

That first step is the object of this paper. Just as in fiat space one must go to two loop 
order to see corrections to the derivative terms, so we will compute the fully renormalized 
self-mass-squared at one and two loop orders. In a later work we will employ the result to 
study the behavior of linearized perturbations in the quantum-corrected field equations. 

In Section II we explain the Schwingcr-Keldysh formalism that one must use for cosmolog- 
ical settings in which the state is specified at a finite time rather than in the asymptotic past 
and future. We also derive the Feynman rules. The one loop counterterms are computed in 
Section III. Sections IV- VI treat the various two loop contributions to the self-mass-squared. 
In Section VII we give the fully renormalized results. Implications are discussed in Section 
VIII. Section IX is an appendix devoted to various technical results used throughout. 
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II. FEYNMAN RULES IN THE SCHWINGER-KELDYSH FORMALISM 



The Schwinger-Keldysh formalism is a technique that makes computing expectation val- 
ues almost as simple as the Feynman rules do for computing in-out matrix elements [84, 85]. 
To sketch the derivation, consider a real scalar field, (f){x) whose Lagrangian (not Lagrangian 
density) at time t is L[(f){t)]. Most quantum field theorists are familiar with the functional 
integral expression for the matrix element of an operator between states whose wave 

functional are given at a starting time s and a last time £, 



$*[(/>(£)] e*/>*^[<^W] ^[0(s)] 



(8) 



The T*-ordering symbol in the matrix element indicates that the operator is time- 

ordered, except that any derivatives are taken outside the time-ordering. Wc can use (8) to 
obtain a similar expression for the matrix element of the anti-time-ordered product of some 
operator 02 [0] in the presence of the reversed states. 



T*{ol 

[d(f)] Csb] ^(f>{i)] e-^/>*^['^Wl **[0(s)] . 



(9) 
(10) 



Now note that summing over a complete set of states $ gives a delta functional. 



^$[0_(£)] $*[(/.+(£)] = 40_(£)-0+(£) 



(11) 



Taking the product of (8) and (10), and using (11), we obtain a functional integral expression 
for the expectation value of any anti-time-ordered operator O2 multiplied by any time- 
ordered operator Oi, 



T*(O2[0])t*(Oi[0]) 



[#_] (5 (£)-(/)+(£) 



X 



(12) 



This is the fundamental relation between the canonical operator formalism and the functional 
integral formalism in the Schwinger-Keldysh formalism. 

The Feynman rules follow from (12) in close analogy to those for in-out matrix elements. 
Because the same field is represented by two different dummy functional variables, (f)±{x), 
the endpoints of lines carry a ± polarity. External lines associated with the operator O2 [0] 
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have - polarity whereas those associated with the operator Oi [(f)] have + polarity. Interaction 
vertices are either all + or all -. Vertices with + polarity are the same as in the usual Feynman 
rules whereas vertices with the - polarity have an additional minus sign. Propagators can 

be ++, -+, +- and --. 

The four propagators can be read off from the fundamental relation (12) when the free 
Lagrangian is substituted for the full one. It is useful to denote canonical expectation values 
in the free theory with a subscript 0. With this convention we see that the ++ propagator is 
just the ordinary Feynman propagator, 



iA++(x; x') - (O|r(0(x)0(x')) = iA{x; x') . 
The other cases are simple to read off and to relate to the Feynman propagator, 
iA_+(x;x')=(n 0(x)(/)(x') n) = e(t-t')iA{x;x')+e{t'-t)\iA(x;x') 



iA+_(x; X ) = {p 
iA__{x;x')^(n 



(f){x')(f){x] 



n 



e{t-t') iA{x]x) +e{t'-t)iA{x;x'), 



T{cl>{x)cl>{x'))\Q)^^[iA{x;x') 



(13) 

(14) 
(15) 
(16) 



Therefore we can get the four propagators of the Schwinger-Keldysh formalism from the 
Feynman propagator once that is known. 

Quantum field theory has much the same structure in the Schwinger-Keldysh formalism 
as for the familiar in-out formalism. In particular, one-particle-irreducible (IPI) diagrams 
mean the same thing and play the same fundamental role. Because each external line can 
take either polarity there are 2^ variations for each A'^-point function of the in-out formalism. 
They generate an effective action that depends upon background fields (f)+{x) and 0-(a;). 
For example, there are four Schwinger-Keldysh 2-point IPI functions, M'l^{x;x'). The ++ 
one is the same as the in-out self-mass-squared and the others are related as the propagators, 

-iM'^_{x;x') = (^-iMl^{x;x')y , -iM^^{x; x') = (^-iMl_{x; x'))* . (17) 

The various self-mass-squared 's enter the effective action as follows [85], 



r[0+,0_] = 5[0+]-5[0_] 

/ ax \ ax 



0( 



2^ i l+0_(a;)M2^(x;x')0+(x') + 0-(x)M!_(x;x')0-(x'; , 

Note that since counterterms derive from the classical actions — S\(^j^ and —S\(^. 
can be no primitive divergences in IPI diagrams of mixed polarity. 



(18) 
there 



6 



The effective field equations of the Schwinger-Keldysh formahsm are obtained by varying 
with respect to either background field and then setting the two fields equal [85] , 



6(p+{x) 



6<j){x) 



J d''x'{Ml^{x; x') + Ml_{x; x')}(f){x') + = . (19) 



Just as in-out matrix elements of the field obey the in-out effective field equations, so 
expectation values of the field obey the Schwinger-Keldysh effective field equations. Unlike 
the in-out effective field equations, the Schwinger-Keldysh effective field equations are causal 
and purely real. That is, although the equations are generally not local, there is no net 
dependence at point x'^ from fields at any point x'^ which lies outside the past light-cone of 
x'^. Hence the initial value problem is well defined. Also, if the initial value data are real, 
the solution remains real. 

Up to this point the discussion has been general. We begin the process of deriving the 
Feynman rules of our model by expressing its Lagrangian (2) in terms of the unrenormalized 
field (p and bare parameters ml and Aq, 

C = -la^^a.^^'^^v^ - ^ml<j>'^ - ^Ao</.^V^ . (20) 

Specializing to the (D-dimensional) locally de Sitter geometry in conformal coordinates (1), 
and introducing the renormalized field cp = 4>/y/Z gives, 

C = ~Zd^<pd,<prj^'^a''-^ - ^Zml<p\'' - ^Z'Ao0"a^ . (21) 

We now make the usual decomposition of bare parameters into renormalized parameters and 
counterterms, 

Z = 1 + 5Z , Zml = + , Z^Xo = X + 5X. (22) 

Note that we define the renormalized mass to be zero. Of course this can only be enforced 
initially; the point of the calculation is to see what happens later under the impact of 
infiationary particle production! 

The final expression of the Lagrangian distinguishes the free part (1st term), the inter- 
action (2nd term) and the counter- interactions (last 3 terms), 

C = -ld^^d,iprj'"'a^-^ - ^Xcp^a^ - \5Zd^^pd,^pr]^"'a^-^ - Urn'^p'^a^ - ^5X(t>^a^ . (23) 
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The propagator derives from the free part and obeys the equation, 

Propagator =^ rj^''d^(a^-^djA{x;x')) = i6^{x - x') . (24) 

The position-space interaction vertex depends upon the four points Xi_4^, 

Interaction =^ —iXa^6^{xi — X2)S^{x2 — X3)S^{xs—X4) . (25) 

The field strength renormahzation depends upon two points x'^ and x''^, 

Field Strength Renormahzation =^ idZT]>"'d^(a^-^dJ^{x-x')) . (26) 

The mass renormahzation also depends upon x^ and x'^, 

Mass Renormahzation =^ —iSm^a^ 5^ {x — x') . (27) 

Of course the coupling constant renormahzation is the same as (25) with A replaced by 5\ 

Coupling Constant Renorm. =^ —i5\a^5^{xi — X2)5^{x2—X3)5^{x3 — Xi) . (28) 

It remains to solve (24) for the propagator. One might think to look for a function of the 
de Sitter invariant geodesic length £(x; x') but there is no such solution. The reason is that 
the free theory possesses no normalizable, de Sitter invariant wave functions, and hence no 
de Sitter invariant states [86]. By allowing iA(x;a;') to depend additionally upon the scale 
factors a = 0(77) and a' = a{r)') one can find a solution which has the homogeneity and 
isotropy of cosmology [2]. It is expressed in terms of a length function we call y++{x; x'), 

y++{x;x') = aa'H^ ||^- (|?7-?7'| = aa'H^Axl^{x;x') , (29) 

When the infinitesimal parameter S vanishes this has a simple relation to the geodesic length, 
y — 4sin^(i7£/2). One can see from (13-16) that the four propagators must all depend the 
same way on a, a' and length functions y±± defined as follows, 

y+_{x;x') = aa'H"^ \^\x-^f-{r]-r]'+i5f\ = aa' Axl_{x; x') , (30) 
y-+(x;x') = aa'H^ [\\x-x'\\^-(r]-r]' -i5)^] = aa'H^Axl^(x;x') , (31) 
y..{x;x') = aa'H^ \^\x-^f -{\r}-i\+i5f\ = aa'H^Axt_{x;x') . (32) 

Because of this close relationship between the four propagators we shall often drop the ± 
subscripts and treat expressions for general y — y{x; x') and Ax'^ — Ax'^{x\ x'). 
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To facilitate dimensional regularization we express the dimension of spacetime in terms 
of its deviation from four, D — A — e. It is also convenient to break the propagator up into 
three terms, 

iA{x; x') = A{x; x') + B{x; x') + C{x; x') , (33) 
A(x-x') = r(l-f) (aaQ-^^i 



r./ A H^'" f 2r(3-f)/|/xf /7rexr(3-e) r(3-e) , , I 



^2— e oo 



1 T{3-e+n) ^yy _ 1 r(3-f +n) n+f 



(36) 



_nr(2-|+n)M^ n+| r(2+n) V4. 

Powers of A{x;x') are the ultimate source of ultraviolet divergences. The term B{x;x') is 
not ultraviolet divergent by itself, 

hmB{x;x') = {^f {-2ln{H'Axy4) - l} = -^^ln{^H'Ax'). (37) 

However, we must retain e 7^ whenever B{x; x') multiplies two or more powers of A{x; x'). 
The term C{x;x') actually vanishes for e = 0. We only need it when multiplying by more 
than two powers of A{x; x'), which can not happen for the maximum number of propagators 
— {iA{x; x')y — we shall encounter in this computation. 

III. ONE LOOP COUNTERTERMS 

Neither the non-fiat geometry, nor the Schwinger-Keldysh formalism, nor even the use 
of position space rather than momentum space, alters the basic structure of perturbative 
quantum field theory. In particular, the diagrams all look the same as in flat momentum 
space, and working at two loops requires knowledge of the one loop counterterms. The 
purpose of this section is to compute the order A contribution to Sm^ and the order 
contribution to SX. As in fiat space there is no one loop contribution to 5Z. 

Fig. 1 depicts the one loop contribution to the scalar self-mass-squared. We can extract 
6m^ at order A by focusing on the ++ term, 

"A 



iM^.Ax;x') 



—iA_^._^_{x; x) + 5w? 
2 



a^5^(a;-a;') + 0(A2) . (38) 



Although the coincidence limit of the massless propagator vanishes when using dimensional 
regularization in flat space, the de Sitter curvature of our background gives a nonzero result 
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FIG. 1: Generic one loop self-mass-squared. 



that can be read off from (33), 

iA{x';x') = \im^tA{x;x') = (^^^§^{2 ln(a) + tt cot(^) | . (39) 



Because of the finite, time-dependent term we cannot make the one loop self-mass-squared 
vanish for all time. Our renormalization condition is that it should be zero at t — 0, which 
implies, 

\H^-' r(3-e) 



Sm — 



7rcotQ+0(A^). 



25-7r2-f r(2-|)" """^2 
The renormalized ++ self-mass-squared is therefore, 

MlJx; x') = ^ ln(a)S\x - x') + O(X^) . 

The — case is minus the complex conjugate. 



(40) 



(41) 



Mi_{x]x') 



87r2 



a''ln{a)6\x-x') + 0{X' 



(42) 



The +- and -+ cases vanish at this order because there are no mixed interactions, 

Ml_(x;x') = -{Ml^{x;x')y ^0 + O{X^) . (43) 








FIG. 2: Generic tree and one loop vertex function. 

We can extract the order contribution to 5 A from the one loop corrections to the 
++++ 4-point IPI function. Fig. 2 shows the diagram topology. A simple application of the 
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Feynman rules gives, 



—iV^{xi, X2, X3, X4) = —i{X+6X)af6^{xi—X2)S^{xi—X3)6^{xi—X4) — 2^^^ 



„DD 



1 2 r 1 2 

iA++(xi; X4) 5'^{xi—X2)S'^{x3—X4)+a^af iA++{xi;x3) S^{xi—X4)S^{x2—X3) 



tA^4xi;x2)] S''{xi-X3)6''{x4-X2)\ + 0{X^) . (44) 



In Appendix B we show that the square of a propagator can be expressed as a local term 
which diverges for e — > 0, plus a nonlocal term which is integrable for e = 0. For the ++ 

propagator the result is equation (95) which we reproduce here, 



(iA++(x;x')) 



2 ' r(l-|) __4^26ri3/_. [aa') fln{fi^Axl_ 



2%2- 



2) ^-4f2exi3 



S'^ix-x') - 



\ Ax 



2 

++ 



+ 



2%4 aa'Ax^^ ^ 2%^ V 4 ^■"++; • ^^^^ 
To fix we are interested only in the local, divergent term. There is such a contribution from 
the s-channel — {i/S.++{xi\Xi)^ — the t-channel — {i/S.++{x-i]X3)^ — and the li-channel 
— {iA.++{xi] X2)) — graphs. They each make the same contribution, so the divergent part 
of the vertex is, 

— il+(xi, X2, ^3, Xi) = + SX)ai~'^5^{xi — X2)5^{xi — X3)5^{xi — X4) 
i3AVr(l-f) 



2%2-f (i-e)e 



af5^{xi — X2)S^{xi — X3)5^{xi — X4) 
X (Nonlocal Finite) + O(A^) . 



(46) 



The simplest renormalization condition is to absorb the entire local, constant part at order 

1677^-2 e(l - e) ^ ^ 

Note that while the divergence is proportional to the fourth power of the scale factor, 

the counterterm is proportional to the Z^-th power. This means that expanding af~^ = 

af{l — eln(ai) + 0(e^)), gives an additional, finite local term. The full renormalized ++++ 

IPI 4-point vertex is, 

/ 3;^2 \ ^2 r 

-iVR{xi, X2, X3, x^)^ -i (a - ^^ln(ai) j a^(5''(xi-X2)5^(xi-X3)(5^(xi-X4) + ^^< 



(0104) di 



ln(/i^Aa; 



Ax^ 



14 



+ 4i7^(aia4)' 



^\^2^ '"^ - H\a,a,)Hn'{y^H'Axl) 



x5^{xi-X2)5'^{x3-X4) + A^(i-channel) + A^ (li-channel) + O(A^) . (48) 
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IV. THE SUNSET DIAGRAM 



The first two loop contribution to the scalar self-mass-squared is known as the sunset (ss) 
diagram. Fig. 3 depicts its topology. A simple application of the Feynman rules gives, 



^ML± = ^(T^A)(T^A)(aa')^-^(^A±±(x;a;'))' . 



(49) 



FIG. 3: Generic Sunset Diagram. 

In Appendix C we show that the cube of a propagator can be expressed as a local term 
which diverges for e — > 0, plus a nonlocal term which is integrable for e = 0. For the ++ 
propagator the result is equation (108) which we reproduce here, 

-'T^(l - f ) {aa')-'^+l'd^5^{x - x') , i3H^~'i2~'a-^+^' ^ 



{iA++{x;x')f ^ 

(i-f)(i-!r(i- 



29 



TT 



4-e 



;i-H(i-e)(i-!e)6 



+ 



(1 - ^e]e 



|) .Hay Tjl-e) 



[aa 



ln{n'^Ax 



2 ^ 



^2^1 ' e 
3H^(aa')-^ 



TTCot (^y ) + 21n(a) 



Axl^ 



+ 



2^6 



In 



Hei X ln(/i2Aa;2 



+ 



In'i^H^Ax^ 



2/i 

2 

++ 



26-^7^4-6 

U^(X-X') 



llnV^ 

^4 Axl+ 



In^^H^Axl^). (50) 



2%6 aa' Axl+ 2^7r^ v 4 

When this is multiplied by (aa')^~^ a few simple rearrangements allow us to recognize the 
spacetime dependent part of the the first term as essentially a field strength renormalization. 



'''+'^d^S''{x-x') = a^+ta2(ai+t5^(x-x')) 

= a^+' - {2+e)Hado - (l + ^) (2 + ^)i^V| S^'ix-x') . (51) 

Doing the same thing for the second term reveals what is essentially a mass renormalization 
plus some finite local terms and an overlapping divergence, 
i3/j-V / 2tt y /(l-|)(l-f )r'(l-f ) / jJqy r(l-e) 



2H^ ^Hi^J \ e2(l-|e) ^ 2^1 ^ 



7rcot(^)-F21n (a) 
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\HuJ I e2(l-|e) 
1. , , 1 



267r4 



e2(l- 



— ln(a) + -ln^(a) + 



3 , /i/ 
--7 + In — 

4 ^ 12/. 



TT cot - 

V 2 

ln(a)U^(a;-a;') . (52) 



In displaying the final result we exhibit the field strength and mass renormalizations, as 
well as the overlapping divergence and the finite, local terms. Our result for the ++ case is, 



M 



r2(i-|) 



SS+ + 



2i037r4V'^ (l-fe)(l-|e)(l-e)e 
21037,4 a^l^\n{a)d^-{2\n{a) + l)Hado-{2\n{a) + ^)H^a^^6\x-x') 
X'H^ ^ 27r V f (l-i)(l-i)r^(l-i) ^ff r(l-6 



- 2HadQ - 2H^a^] 5^ {x-x') 

3^ 



-f— 

\Hn) 



277r4 Vii-^. 



6^(1-16) 



A2i/2 27r X e ln(a) 4 .^^^ , ^ 1 



iX^aa' 



24 



In (//^Ax^^) 
. Ax2^ 



2^4 I 2 

-H^aa'd'^ 



2/1^ e 

^ ln^(a) + 



-TT 



cot(yUa^5^(a;-a;') 



ln(a) 



^ 2/1 ) 



^xl. 



4Ax2^ 



(63) 



The — case is minus the complex conjugate. 



r2(i- 



-2//a9o -2/^2021 (^^(x-x') 



A^ 



2103774 



|ln(a)a2 - (2 ln(a) + l)HadQ - (2 ln(a) + ^H^a'\5\x-x') 



X'H' ^2nyi (1-f )(l-f )r2(l-|) ^ r(l-6) 



-TT cot 



+ 



iA2aa' f 1 



2%^ 1 24 



In (^fi'^Axl 
Ax2_ 

-HHao! 



2^77^ I 2 



ln(a) la^5^(a;-x') 



Hel ^ In (//2Ax!_) In' (//'Ax2__) 
A^2 + 4Ax2 



(54) 



The mixed polarities lack the local terms and differ by a sign, 



[aa 



[iA+_(a;;a;')] 



/m3 



A2aa' f 1 



2%6 24 



-H^aa'd^ 



In 



'i7et\ ln(//2A4_) 1 In' (/i2Ax 
A4I ^4 A^T 



In (^iJ,'^Ax'^_ 



Axl. 



H\aa') 



y{^H-Axl_) 



Ax\_ 
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+:_(aa')Hn^(^H^Axl_]\ , (55) 



In (//^Aj;?_+) 



2/i 



Ax? 



■+ 



Aa;2 



H\aa') 



+^(aa')' In^ f^i/^^xi/ 



(56) 



V. THE SNOWMAN DIAGRAM 



The primitive two loop graph is known as the snowman (sm) diagram. Unhke the sunset 
diagram, it contributes only to M^_^ and M^_. Its generic topology is depicted in Fig. 4. 
The right hand diagram depicts the 1-loop mass counterterm which naturally combines with 
it. 




FIG. 4: Generic Snowman Diagram with mass counterterm. 



Because both external points of this diagram must have the same polarity we denote the 
contributions with a single subscript sign, Af^^ and M^^^ . The internal vertices are summed 
over both polarities, as always in the Schwinger-Keldysh formalism. A simple application 
of the Feynman rules gives. 



-iMi 



D 



5^ (a; - x')l (iVa"^<' [iA++{x;x")f - [iA+_(x; x")]' 



X ■ 



iA{x"-x")-i5m'^\ . (57) 



Since both diagrams of Fig. 4 have the same lower loop, they possess a common factor in 
the first curly bracket. The first term within the final curly bracket derives from the left 
hand diagram, while the second term comes from the right hand diagram. 
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The coincident propagator and the mass counter term are calculated in (39) and (40) 
respectively. Once they are inserted into the above equation, one finds 



-iM 



2 

sm-f 



\^H^ ^a^ ^r(3 e).5i?(^_^YdVa''Mn(a'')|[iA++(x;a:'')]'-[iA+_(x;x'')]' 



25-^7r2-f r(2-|) 



(58) 



The squares of the various propagators are computed in Appendix B and the results we need 

are (95) and (96) 

iA2i7V/27rx.r(l-e). , 



-iM 



2 

sm-|- 



■f— )^ 
\Hn.J 



e 



■ln{a)S^{x-x') 



210^6 



ln(/x2Aa;^_^) \n{^^^x'^ 



Aa;^_^ Ax^ 



2%6 



6\x-x'){-dl) j dVa"Mn(a") 

5\x-x')jd^x"a"^Ha")^ ' ^ ' ^ 



Ax2 



Aa;2 



-^5\x-x')Jd^x"a"nn{a") 



210 

g^g^ / 27r y / ln( 



ln<^A.^Vln<^A.^ 



(59) 



2%^ 

2107^6 



f7ln(a)U^(^^(a;-x') 



+ 



S\x-x'){li + AH^ah - H^a^h) . 



(60) 



The integrals Ji, I2 and Is arc evaluated in Appendix D and the results are expressions 
(119), (125), and (129). Substituting these and adding everything up gives. 



■— r— r -rr- (x — x)-^ =— - a (x—x < m (a) 



+ ^ln'(a) + 



\ ^H\ S 2 21 w ^ 238 13 2 4^,^, 

n( — H 7--7r^ n(a) \ 7r^ + -C(3) + 

^2/1^ 3 ' 9 i ^ ' 81 54 3^^ ^ 



81 



n=l + 1)' n=l ' ^ 



(61) 



(n + 2)3 ^^n(n + 3)3j 

Note that the overlapping divergence from the first term is the same as the one from the 
sunset diagram. The — case is minus (61). 



VI. THE COUNTER DIAGRAMS 



We have already encountered the mass counterterm in Section 111. At the order we are 
working there is also a one loop diagram with A replaced by SX, and there is a field strength 
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5Z 

— e — 

FIG. 5: Generic Counter Diagrams. The one on the left is the same as Fig. 1 but with SX in place 
of A. The one on the right is the field strength renormalization. 

renormalization. The generic diagram topologies are depicted in Fig. 5. Note that neither 
diagram contributes to the mixed polarities. 

The contribution of the one loop counter vertex diagram to the ++ polarity is, 

-iMl^ = -i^a''iA{x';x')5''{x-x') . (62) 

Of course we can read off the answer from Section III, 

Ml, = i.-^ ^^|21n(a) +7rcot(-)|^^(a;-a;0 , (63) 

i72A2 . 27rNe f3r(l-e) /vrex 3,2/xl x.n, 
~ 2%J-(^) |2-4^''<=°*(y) --^iAa)]a*S"(.-.') . (64) 

The — case is minus this. 

The contribution of field strength renormalization to the ++ polarity is, 

Mf^, = -SZr]'"'d^(a^-^dJ^{x-x')) = -5Za^-'(d'' -{2-e)Had^5^{x-x') . (65) 

The only primitive graph with a divergence proportional to a^(d'^—2Hado^5^{x — x') is the 
sunset diagram. Comparison with our result (53) for it implies, 

(^Y r^(i-t) .... 

2^W^iJ?^ (l-fe)(l-|e)(l-e)e ■ ^ ^ 

Although some finite, local terms result from expanding — 1 — eln(a) + O(e^) there are 
no overlapping divergences from the field strength renormalization at two loop order. 

Ml, = -^-(^Y- -aHd^-2Hado)5''{x-x') 

iz+ 2^037r'iV^^ (l-fe)(l-ie)(l-e)e ^ ^) \ ) 

+^I^a2 [-ln(a)a2+(21n(a) + l)Mo] 5\x-x') . (67) 
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After field strength renormalization, all remaining divergences reside on a^5-^(x— x'). Such 
contributions comes from the sunset diagram (53) and the counter vertex (64). Comparison 
with these results suggests that the simplest choice for the two loop mass counterterm is, 



A2 



A^g^ / 27ry i (Hy 

rri-e 



+ : 



H 



12(l-|6)(l-fe)(l-e)e ' V2^) 



:i-i)(i-fr(i-i: 



+ 



2e 



/7re\ i 

-TTCOt^yj + 



1 

16 



238 13 2 4^, , 



(68) 



Of course the mass counterterm is Mf^_^ — 5m?a^ ^5^{x — x'). Using the expansion a ^ — 
1 — eln(a) + |e^W(a)+0(e^) we see that the two loop contribution is, 



A2 



2e 



277r4 

-TTCOt [— 



r^fi-^ 



(i_l)(i-|)r2(i_l) 



l-|e)e2 



277r4 



In'(a)- 



2 y 

5 



V2/i^ 12(l-|e)(l-|e)(l-6)e ^2^1^ ^ 



r(l-e) ,f7re^\ 4-^, . , H\ 27r y ln{a) 



ln(a) + 



1 



238 13 2 4^ 1 4^4( _ (gg^ 

16 81 54 3^^ M ^ V ^ 



Note that the nonlocal divergence precisely cancels the nonlocal divergences from the sunset 
diagram (53) and the snowman diagram (61). As usual, the — polarity is minus (69). 

VII. THE FULLY RENORMALIZED RESULT 

To get the full two loop result for the ++ case we add the contributions from the sunset 
diagram (53), the snowman diagram (61), the counter vertex diagram (64), the field strength 
renormalization (67) and the two loop mass renormalization (69), 



M 



2+ 

i\ I aa 



' ';m4- ' r\r-\- ' A /, A- ' f\'rr\ 



2%6 I 24 



In (^/i^Ax^ 



'(5Z+ 

r2/„„/\2o2 



^^~2^' 



+ 



In^ (/i^Ax- 



4Aa;2^ 



-H\aa') 



Axl^ 



+ 



293 TT^ 

X'H^ { 4 



6 ^ ' V 4 
ln(a) + (2 ln(a) + 1) Hado ^5^{x-x') 
23 



, , ,--ln3(a) ln2(a) + 

2^77^ 9 ^ ^ 18 ^ ^ 



13 ^ / H \ 2 o 
— +31n( — )--7r2 
3 ^2fjL^ 9 

(n+2) 



ln(a)|a^(5^(x-x') 

■(n+3) 



X2H"2 f „-3 00 I r 00 „-(n+2 00 „-(n+3) ^ 
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As with the renormahzed stress-energy tensor [3] , we conjecture that the exponentially falling 
terms on the last line could be absorbed into a perturbative modification of the initial state. 
The — polarity is minus the complex conjugate, 



iX^ { aa' 



2%6 24 



In (^jji^Axl 



-H'\aa!fd'' 



i/etx In (/x'Aa;!_) In^ (/x^Ax? 



In(^) 



4Aa;2 



A' 



-J a'l - ln(a)9'+ (2 ln{a) + l)Hado \S\x-x') 



293 

X'^H^ r 4,3,, 23, 2, ^ 



13 



+ 31nf^)-^7r- 



2//^ 9 



In (a) UV(x-x') 



X^H^ f ^ n + 5 
1^ ' — 



oo „-(n+2) oo „-(n+3) ~1 

^'^"^^+4E 7?r-^+4E i:^:^ [aV(x-x') .(7i) 



As usual, the +- and -+ polarities have the opposite sign and lack the local terms. 



2+- 



2^6 24 



In (/x2Aa;2_) 



Ax2 



-H'(aa'?d^ 



Hek In (/^'Aa;2_) In^ (/x^Aa;^ _) 



^ 2^ ^ 



Ax"" 



4Ax2 



Aa;2 



(72) 



2%6 1 24 



In [jj^Axl 
A^ 



-H\aa'fd^ 



Hek In (/i'Ax!^) In^ (/x^Aa;!^) 



V 2/x >' 



Ax^ 



4Ax2 



(73) 



VIII. DISCUSSION 



We have computed and fully renormalized the four Schwinger-Keldysh self-mass-squared's 
at one and two loop orders for massless, minimally coupled Lp^ (2) on a locally dc Sitter 
background (1). The next step is using our results to study the quantum-corrected, effective 
field equations at linearized order (7). The term "M^(x;x')" called for in (7) is the sum of 
Ml_^{x;x') and Ml_{x;x'). The analysis — which will be presented in a subsequent paper 
— is very similar to that already done for de Sitter photons using the one loop vacuum 
polarization [87] and for de Sitter fermions using the one loop self-energy [88]. 
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There are two interesting issues to investigate. The first is stabihty. Secular effects in 
this system are driven by inflationary particle production which forces the scalar up its </7^ 
potential. One would expect this to generate a positive mass-squared, and it is significant 
that our one loop result (41) is positive. The two loop result (70) is a little too complicated 
to interpret without a detailed calculation of the type proposed, but we expect no problems 
with stability. 

The second interesting issue is an apparent paradox in the stochastic solution of this model 
which was given in 1994 by Starobinskii and Yokoyama [89]. As we noted in the introduction, 
this model is stable for two reasons. First, as infiationary particle production forces the 
scalar up its potential the classical force pushes it back down. Second, the curvature of 
the potential ought to induce a positive mass-squared which should reduce the inflationary 
particle production. However, the solution of Starobinskii and Yokoyama implements only 
the flrst effect. Although it would be possible to implement the second effect as well within 
the context of Starobinskii's stochastic technique [90], including just the flrst effect gives 
perfect agreement at leading logarithm order as high as this has been checked. The question 
is why should this be so? A possible answer is that the nonzero self-mass-squared remains 
always down by one factor of ln(a). Whereas the leading logarithms of the energy density 
give a series in Aln^(a), 




It is worth noting that, had gravity been quantized, there would have been similar cor- 
rections from graviton exchange. While the physics of the pure scalar model pretty much 
guarantees that a positive self-mass-squared must be induced, there is no such argument for 
gravity. If the induced self-mass-squared turns out to be tachyonic it might have important 
implications for scalar-driven inflation. Further, there is no reason the gravitational process 
should ever stop. The graviton should not experience a classical force since it has no non- 
derivative potential. Nor should the graviton ever develop a mass. This might be signiflcant 
for scalar-driven models which inflate for a very large number of e-foldings. 

We close with a comment on accuracy. Two loop computations are generally intricate 
enough that this becomes an issue. In this regard it is worth noting that three distinct 





(74) 



it may be that the self-mass-squared behaves as such a series times A ln(a) 



M\x;x')[^^^^,^^ = H''aH\x-x') x Aln(a) 5^m^(Aln'(a)) 



(75) 
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diagrams conspire to cancel the overlapping divergences. This is significant because there 
are no overlapping divergences at this order in massless c/?^ on fiat space. 

It is also a significant correspondence check that we get the same results for 5X (at 
one loop) and 6Z (at two loops) as for massless ip'^ in flat space. This had to be the 
case because the divergences these counterterms absorb are formally logarithmic. Because 
geometric corrections make divergences more finite, any geometric corrections to 5 A and SZ 
must be finite. The same is not the case for Srn^ because the divergences it cancels are 
formally quadratic. Indeed, Sm^ vanishes for massless (p^ in fiat space whereas we found 
logarithmically divergent corrections proportional to H'^ at both one and two loop orders. 



IX. APPENDIX 



A. Partial Integration Identities 

In almost all cases derivatives of the conformal coordinate intervals Ax'^j^{x;x') can be 
taken generically, without regard to polarity and without worrying about delta functions, 

/ 1 \ (n—me](n — 2 — (m — l)e] , . 

aM— ^ ^ (n^2orm^l). (76) 

The exception is l/Ax^"*^, 



r(i-f) ^ ' ' V^4-V V^^-+V ' 

(77) 



These identities allow one to extract more and more derivatives of inverse powers of Ax 
until the result is integrable in D — A dimensions. For m 7^ 1 we have, 

1 92 / 1 \ 



2 

±± 



A^--^ (^2-me)(-(m-l)e) V^^±± 

1 (I 



2-mej > (78) 
^ ' ^ (79) 



At this stage one adds zero in the form of the appropriate identity (77) and then expands 
the nonlocal term in powers of e. For example, the ++ term gives, 

\^Ax2-"»V yAxl-"^' Axl-' I r(l-f) V y> \ ) 



+ r(i-f) ^ 



The cases we require are, 



A4-2^ (2-2e)(-e) r(l- 



+ 



(2-2e)(-e) \ A^^ 



A4-3^ (2-36)(-26) r(i-i) 



+ 



(2-3e)(-2e) I A^^ 



A4-3^ (4-3e)(2-3e)(2-2e)(-2e) r(l-f; 

^(4-3e)(2-3e)(2-2e)(-2e) \ Aa;2^ 
The — variations differ only in the sign of the delta function term, 



Axli^e (2-2e)(-e) r(l-f) 



+ 



(2-2e)(-e) I /\xl_ 

S^'ix-x') 



Axli^e (2-3e)(-2e) r(l-f) 

^-3e^2 / eln(/x2Aa;2_) + [iJ? /^xl_) + 0(e3 



(2-3e)(-2e) \^ AxL 

1 ^47r2-f 2^ _ 

Ax6-3e (4_3e)(2-3e)(2-2e)(-2e) r(l-f) ^ 

^-36^4 /eln(/x2Aa;!_) +0(e2)^ 

^(4-3e)(2-3e)(2-2e)(-2e) A^^^ 

And the mixed terms have no delta functions, 

1 _ ^\e\n{^i■'i^xl_) + leHn\^J?/S.xl_)+0{e^)' 

A^E^ ~ (2-2e)(-e) I Ax2_ 
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(89) 
(90) 



Ax4-3^ (2-3e)(-2e) \ Axl_ 

Ax6-3^ ~ (4-3e)(2-3e)(2-2e)(-2e) y Axl_ 
Replacing +- with -+ gives the final case. 

B. The Square of a Propagator 

When the square of a propagator multiphes regular expressions we first decompose it 
according to (34) and then take e to zero in all but the term, 

(iA{x;x')y = (^A{x;x')y + 2A{x;x')(^B{x;x') + C{x;x')) + (^B{x-,x')+C{x;x')y, (91) 

_ T^^iaar^ _ H^lni^H^ + J^^n^^H^Ax^) (92) 
2%4-e 2%4 aa'Ax^ ^ K a r ^ ' 

We then employ the appropriate identity — (82), (85) or (88) — and take e = in the 
nonlocal part. The squares of the +- and +- propagators are finite, 

V + J OPiT^ A^2 94„1 nni\^^ ^Stt^ V 4 ^ > 



+ 



2 

+ 



(94) 



V \ Axl^ ) 2%4 aa'Axl^ V 4 

The square of the ++ propagator includes a local divergent term in addition to the corre- 
sponding ++ finite terms. 



K,(.;.'))' ^ -4^^^a-^-5-(x-x')-^a^ fi^^ 
^ ' 2%2-2 (l-e)e ^ ' \ Axl^ 



The square of the — propagator is the complex conjugate, 
UA (x-x')Y ^/i^- r(l-f) ^ (ggQ-^ / ln(^Wj 
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C. The Cube of a Propagator 



When the cube of a propagator multiphes a regular expression we again decompose it 
according to (34) but now we must retain e 7^ in the and SA^B terms, 



(iA{x; x')y ^A^ + sA'iB+C) + 3A{B+Cf + {B+C)\ 

r3(i-f)(aa')-'+i^ 3/i'2-T(i-|)r(2-e) r 2r(3-f)r(2- 



(97) 



2671^-1^ Ax' 



,6-3e 



27- 



+ 



7rc()t(^) + lii(aa') 



+ 



%aa'y 
In 



r(3-e) Ax^-^' 



\ A 



2%6 aa'Ax"^ 2%^ - v 4 

The term is easy to read off from the identities (84), (87) and (90). For example, the 
case gives. 



(98) 



[aa 



n-3 



ln(/^^Aa;;_) 



Axl_ 



(99) 



267^6- Ax^-_^' 2ii7r6 
The ?,A^B term is more complicated on account of the 1/e divergence in its prefactors. One 

must first partially integrate and multiply before taking e = in the nonlocal term. For the 
H — case we have. 



3(^'S)+_ = 



27-^7r6-5^ 



r(i- -) 
^ 2' 



X • 



0(1 



t)r^(l-l) rH^aa\i ^ 
V 4/^2 



hi{ii^Ax\_) ^ ln^(/i2Aa;^. 



Axi 



Axi 



r(i-6) 



TT cot ( + In(aa') 



ln(/i2Ax^_) , 3 ln='(/i2Ax2 



Axt 



Ax"^ 



+- 



At this point we digress to recall some basic series expansions. 



r(l-e) = l + 7e + 



12 



/7re\ 2 

TTCOt = - 

V 2 ^ e 



12 ^ ' 



It follows that the various prefactors in (100) are, 

(i-i)(i-tr(i-i) ^ 

2(l-|e)e 



■— ^ ^TTCOt — 

4 V 2 

r(i-e) 



l + 2 + 3+o(e), 
2e 2 8 ^ ^ ' 

-i + 0(.). 



Making use of these expansions in (100) gives. 



3{A^By 



3H^{aa') 
2%^ 



Hel X \ii{ii''Axl_) 1 Infill'' Ax\_) 



V 2/x ^ 



Axl_ 



+ 



Axi_ 



(100) 
(101) 

(102) 

(103) 
(104) 

(105) 
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Employing (99) and (105) in (98) gives 



+ 



3H^{aa 



J\-2 



+ 4 



lln='(/x2^4_) 



Ax? 



2%6 

2%6 



aa 



The -+ case follows by complex conjugation, 

ln(/x^Aa;^. 



(^A_+(a;;a;'))-^ 



2ii7r6 



Ax2_ 



+ 



3H'^{aa 



l\-2 



2%6 



■ /ifetx ln(u2Axi_) 

/n^) A4_ 



4 Ax?.+ J 2%6 aa'Ax?. + 2971^ V 4 +) ^ ' 

;ator includes a local divergent term in addition to the correspond- 



The cube of the ++ propags 
ing ++ finite terms, 



{iA^4x;x')r 

(l-f)(l-|)r^(l-f) .Hay _ 
(l-|e)e2 V2«>' 



29 71^-' 

2/x 



(1-|6)(1-6)(1-|6)6 



-4+2e 



r(i-6) 



TT 



cot (^y j + 2 In (a) 



26-67^4-6 
5^(x - x') 



2i%6 



ln(//^Aa; 



Ax2^ 



+ 



?,H\aa')-^ ^ 
2%6 



/ifetx ln(/x2Aa;^+) lln^(/i2Ax 
^ Axi.. ^ 4 



Ax\+ 



3H^ ln^{'§H^Ax 



,2 

++ 



2%6 



aa 



Axi+ 

|^ln^(f^^Axg.(108) 



The cube of the — propagator is the complex conjugate, 

i/x-"T^(l - f ) (QQ^)-3+i^a2(^^(x - x') i3H^-'pL-'a-^+^' [ 
29 TT^-^ (l_|e)(l_e)(l_|e)e 26-^77^-^ \ 



iA__(x;a;'))^ 



(l-0(l-T)r^(l-§) /:^y _ 
(1 - |e)e2 V J 



71 cot (^y ^ + 2 In (a) 



2ii7r6 



'ln(//-A,7-!_ 
A^2— 



+ 



2%6 



/Jeix ln(//2A,r2_) ^ 1 hi-(//2A,r 



V 2/x ^ 

2 



Ax? 



^2%6 aa'Ax?._ 



Ax?._ ' 4 



297r6 



D. Evaluating the Integrals for the Snowman Diagram 



The first integral we need in equation (59) is, 

"ln(/i2Ax^+) 



h = -dlj ciVa"2ln(a") 



ln(/i2Ax^ 
Axl_ 



(110) 
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The first step is to make use of the identity, 



In^ {/^x^)-2hi (ij,^Ax') 



(111) 



Ax2 8 

Because the derivative is with respect to we can move it outside the integration over x"^. 
The remaining integrand possesses only logarithmic singularities. 
We define the coordinate separations, 



Af] = Tj — Tj" 



1 1 -> ->// 1 

r = \\x — X 



;ii2) 



The ++ and +- terms cancel for A77 < so we can restrict the integration to A77 > 0. In this 
case the logarithms can be expanded as. 



In 



li^Axl^ 



In 



li\Ar)^ - r^) ± i7r0{Ar)^ - r') , 



(113) 



We then perform the angular and radial integrations and act the first three derivatives. 



rv r^v 
h = i27r2a^ / dr]"a"Hn{a;') 

Jvi Jo 



= i27r^d^ r dr]"a"Hn{a")A7]^ 

Jrii 

= iSTT^ao r dr]"a"^ In {a") In (2/XA77) . 



2 2\ 

r , 



-In (2/XA77)- — 



(114) 



The lower limit rji — —H"^ (t" = 0) derives from the fact that the state is released in free 
Bunch-Davies vacuum at this instant. 

To perform the remaining temporal integration first change variables from rj" to a" — 
— l/{Hr]") and make the corresponding change in the derivative, 

d r , „, . , r2/i/ 1 1^ 



h = ^87^V^ r da" ln{a") In 



;ii5) 



Then break up the second logarithm and act the derivative on the two nonsingular terms. 



ln2(a) - In f ln(a) - T da" \n(a") In f 1 - — 
\H y da Ji \ a 



Now substitute the expansion. 



and integrate termwise. 



/ da"\nia") In 1 - — U -a ln(a) - 2 + — + V 



(116) 



(117) 



(118) 
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Combining terms gives the following final result, 

ln2(a) - In 



2/i 



ln(a)+Ma)-l + --i:|^ 



(119) 



The second integral in (59) is, 



J d'x"a"Hn{a") 



(120) 



Proceeding as before we extract a derivative with identity (111) and then exploit (113) to 
obtain. 



2 J V^v^ 



r^) 



di'a"Hn{a"^ 



Art / 

drr^ In / a „2 _2^ 







(Ar/" - r') 



(121) 



Now perform the angular integrations, make the change of variable r = A77 z, act the 
temporal derivatives and change variables from 77" to a", 

i/2 



72 = -i27r' r dr^^a"^ ln(a") dzz'^d^ 

Jrii Jo 

= -i27r^r dri"a"^ \n(a")Ar] dzz^ \ 6 In 

■/r/i JO [ 

!f.a"ln(a")(l-041n(l 



Ar73ln(^^Ar;^(l-z2^ 



+ 7 



7/2 



A small rearrangement gives. 



1 1 

a 



(122) 



da"{l 31n(a") +41n'(a") - 4/ da" [1 In (a") In 1 



The only non-trivial integral left in (123) can be evaluated using (117) 



(123) 



da"[ 1 - ^ )ln(a")ln(l - ^ ) =-a{ ^Ha)-l + ^+j: 



2 00 



a J \ a ^ 

The final result for I2 is therefore. 



12 



n=l 



-(n+1) 



-(n+2) 



n{n+iy n{n + 2y 



(124) 



2 In^(a) - 1 ln(a) + ^ + ^ - 5a-^ + ^^^^ +4E 



-(n+l) 



-(n+2) 



n(n + 1)2 n(n + 2)^ 



■ (125) 



The initial reduction of the last integral in (59) is familiar from li and 
h = j d'x"a"' \n{a") 



26 



= a67r2 /' d^z' r cir/"a"^ ln(a") Ar/^ In 
JO J\ a" \ a , 



il67r2 r^da" 



In 



13 
18 



(126) 



At this stage it is best to integrate by parts on the term ln(a")/a" = d/da" {\v?{a")/'2}j^ 



1 /13 



a V12 



Ma")) 



The first term is straightforward and the final term can be evaluated using (117), 

0- 9' '"V). (: - g ^ - i . . |C(3) . g 4) in,) 

o r 9^-('"+2) 1 1 

■ (128) 



-1 OO 

--ln» + 2E 

y n=l 



n(n + 1)^ n{n + 2)^ n(n + 3)^ 
Substituting (128) in (127) and adding everything up gives, 
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